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^ ■ ABSTRACT 

in '. 

t-H ■ Context. Recent 3-D analysis of the solar spectrum data suggests a significant change of the solar chemical composition. This may affect the 
' temporal evolution of the surface abundance of light elements since the extension of the convective envelope is largely affected by the internal 
^^ , opacity value. 

>0 ■ Aims. We analyse the influence of the adopted solar mixture on the opacity in the convective envelope of pre-main sequence (PMS) stars and 
thus on PMS lithium depletion. The surface Li abundance depends on the relative efficiency of several processes, some of them still not known 
r~| with the required precision; this paper thus analyses one of the aspects of this "puzzle". 

O | Methods. Focusing on PMS evolution, where the largest amount of Li burning occurs, we computed stellar models for three selected masses 

JL (0.8, 1.0 and 1.2 Mq, with Z = 0.013, Y=0.27, ff=1.9) by varying the chemical mixture, that is the internal element distribution in Z. We 

analysed the contribution of the single elements to the opacity at the temperatures and densities of interest for Li depletion. Several mixtures 



were obtained by varying the abundance of the most important elements one at a time; we then calculated the corresponding PMS Li abundance 

evolution. 

Results. We found that a mixture variation does change the Li abundance: at fixed total metallicity, the Li depletion increases when increasing 

the fraction of elements heavier than O. 
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1 . Introduction and for stars in open clusters of different ages (see, e.g. Jeffries 

I2000 and references therein). 
In the last twenty years several observational studies of lithium 

( 7 Li) in stars have been performed, many of which focus on The task is not eas y because the tem P oral evolution ° f 
open cluster objects in the zero-age main sequence (ZAMS) the surface chemical abundance (when the accretion phase 
and in the main sequence (MS) phase (see the reviews by has ended ) de P ends on the relative efficiency of several pro- 
Jeffries |2000| Randich M and the references quoted in cesses: mixin 8' microscopic diffusion, radiative levitation and 
Table 1 of Sestito & Randich HUU3). The temporal evolu- mass loss - In Particular, the abundances of the light elements, 
tion of the surface Li abundance has also been investigated which are easil y burned ' are extremely sensitive to the ex- 
theoretically (e.g. D'Antona & Mazzitelli QSZZI review by tent of envelo P e mixin §- There are two main sources of un " 
Pinsonneault 1997 and references therein; Siess, Dufour & certainties in these processes: the first one related to macro- 
Forestini 2000.. The interest in Li depletion in stars comes sco P ic P^cesses and the other to microscopic physical ingre- 
from the fact that, due to its relatively low burning temperature dients ado P ted ' Regarding the former, it is well known that due 
(~ 2.5X10 6 K), Li is a key tracer of the efficiency of mixing to the lack of a self-consistent treatment of convection, stellar 
processes in stellar envelopes. modelers are compelled to use relatively crude approximations 

T .. r. , ff , ,, ,, . ,.„ , , . depending on free parameters, such as the mixing length for- 

In spite of great efforts, the theory is still unable to repro- * . & & 

j ., , iuu- e e t-uj a ■ malism. Even a precise physical treatment of the efficiency of 

duce the temporal behavior of surface Li abundances during F F J J 

., • c /ni\/ro\ j iv/rc u • a ^ the mass loss is not available. The situation is even worse if 
the pre-main Sequence (PMS) and MS phases, i.e. model pre- 

j. .• , . , ., , , ■ e ., c the assumption of spherical symmetry is relaxed allowing for 
dictions do not match the measured abundances for the Sun F F J J & 
rotation, or if magnetic fields are present. Once a scheme has 

Send offprint requests to: P. Sestito, email:sestito@arcetri. astro.it been chosen for the treatment of a given macroscopic process 
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(i.e. convective mixing), the prediction of its efficiency is still 
affected by uncertainties in the basic physical inputs such as 
the nuclear reaction rates, the equation of state and the opacity 
of the stellar matter (see e.g., D'Antona & Mazzitelli 1997 



Pinsonneault 1997. Turcotte & Christensen-Dalsgaard 1998 
Palla l2000l Jeffries 2006 for some results of evolutionary mod- 
els). The canonical scenario assumes spherical symmetry with 
constant mass and adopts the classical Schwarzschild criterion 
for the identification of convective boundaries. 

A fundamental role in Li evolution is played by uncertain- 
ties in the radiative opacity (which affects the depth of the con- 
vective envelope) and thus by the metal content, both the global 
abundance (i.e. the metallicity) and the distribution of the ele- 
ments in the mixture. Recent analysis of solar spectroscopic 
data using 3 -dimensional hydrodynamic atmospheric models 
(see Asplund et al. 2004 Asplund, Grevesse, & Sauval 2005 
- hereafter As04, As05 - and references therein) have reduced 
the derived abundances of CNO and other heavy elements with 
respect to previous estimates (Grevesse & Sauval 1998, here- 
after GS98). GS98 improved the mixture by Grevesse & Noels 
(1993, hereafter GN93), widely adopted in the literature, re- 
vising the CNO and Ne abundance and confirming the very 
good agreement between the new photospheric and meteoric 
results for iron. The Z/X solar value decreased from the GN93 
value (Z/X) =O.O245 to the GS98 estimate (Z/X) =O.O23O and 
then to the values (Z/X) =O.O176, 0.0165 found by As04 and 
As05, respectively. The corresponding solar metallicity values 
(i.e. assuming [Fe/H]=0, see also Sect. [3 for the different mix- 
tures are Z -0.018 (GN93 and GS98), Z -0.013 (As04) and 
Z ~0.012 (As05). However, the problem of a precise determi- 
nation of the solar mixture is still affected by theoretical and 
observational uncertainties and the quoted error on the Z/X so- 
lar value is still about 10% (see e.g. Bahcall & Serenelli 2005). 

In this paper we analyse the influence of the adopted solar 
mixture on the opacity at the base of the convective envelope of 
pre-main sequence stars and thus on Li abundance in the con- 
vective envelope. We focus on PMS since, according to classi- 
cal models, the largest amount of surface Li destruction for low 
mass stars occurs during this stage. The results are understood 
in terms of the contribution of each element to the total opacity 
at the physical conditions typical of the bottom of the convec- 
tive zone during PMS. Clearly this is not an attempt to solve the 
"Li problem", which is much more complex, but a quantitative 
investigation of one of the aspects of the "puzzle". 

A few studies have already been carried out concerning 
the dependence of Li depletion on element abundances and/or 
on the surface stellar opacities for solar-type stars (see e.g. 
Swenson et al. 1994a 1994b and references therein; Piau & 



Turck-Chieze 2002 and Piau, Randich & Palla 2003 1. Swenson 
et al. explored the effect of the opacity values (and of other 
physical inputs) on PMS Li depletion, as well as the effects of 
a variation of the oxygen abundance, trying to match the Li 
abundances observed for the Sun and the Hyades (under the 
assumption of no MS depletion), finding that an enhancement 
in [O/Fe] leads to a larger Li destruction. They also noted that 
heavier species such as Fe, Ne, Mg and Si might significantly 
contribute to the opacity. 



Similarly, Piau & Turck-Chieze (2002i presented PMS 
models with Pleiades (solar) and Hyades-like compositions. 
They investigated the dependence of Li depletion on metal and 
He fractions, by varying the global mixture (and changing the 
total Z); they showed that Li abundances are very sensitive to 
composition, suggesting that the dispersion in the observed val- 
ues in open clusters could be related to small changes in the 
metal mixture among the cluster stars. 

Piau et al. (2003) investigated the sensitivity of the exten- 
sion of the convective zone during the MS for a mass of ~ 1 
M , with respect to the surface content of CNO. Their models 
are calculated by changing [CNO/Fe] while maintaining [Fe/H] 
and the effective temperature unaltered. Thus, the total Z in- 
creases if CNO increases (and vice versa). Under these con- 
ditions, Piau et al. found that an enhancement of CNO results 
into a smaller amount of Li depletion. Piau ( 2005 1 investigated 
the history of Li isotopes in Population n dwarfs, finding a cor- 
relation between the scatter in Li abundance and variations in 
[Fe/O] ratios, i.e. changes in [Fe/O] affect Li depletion more 
than a change in the global metallicity. 

The paper is organized as follows: Sect. 2 briefly describes 
the evolutionary code and the physical inputs adopted for the 
calculations; in Sect. 3 we present the results of our computa- 
tions for models with different mixtures and fixed metallicity 
(Z = 0.013), which are interpreted in Sect. 4. In Sect. 5 we also 
discuss the effect of changing both the mixture and the total 
metallicity value Z. A summary (Sect. 6) closes the paper. 

2. Evolutionary code and input physics 

The models have been computed with an updated version of 
the FRANEC evolutionary code (see e.g. Chieffi & Straniero 
1989 1; the physical inputs adopted for the calculations have 
been listed in the Appendix of Cariulo, DeglTnnocenti & 
Castellani (2004) - see also Ciacio, DeglTnnocenti & Ricci 
(1997). With respect to Cariulo et al. (2004) we updated the 
equation of state (OPAL EOS_2001 1 , see also Rogers (2001), 
Iglesias & Rogers ( 1996 1 and the conductive opacity (Potekhin 
1999 ). The nuclear cross sections for Li burning are taken from 
the NACRE compilation (Angulo et al. 1999). The adopted ini- 
tial value for Li in disk stars is log «(Li)=3.3 2 (see, e.g. Jeffries 

2006),orJV Li/ W H ~ 2x 10 ~ 9 - 

The opacity values are fundamental for model computa- 
tion; we adopt the radiative opacities by the Livermore group 
(Iglesias & Rogers 1996 1 for temperatures higher than 12000 K 
and those by Alexander & Ferguson ( 1994) for the lower tem- 
peratures. As it will be discussed in the next section, we used 
several different mixtures for the high temperature opacities, 
while the Grevesse & Anders ( 1991 ) mixture has been adopted 
for the molecular atmospheric opacities. We also notice that the 
adopted EOS is calculated for the Grevesse & Anders (1991) 
mixture; to our knowledge EOS tables with an updated mixture 
are not available, however, as stated in the OPAL EOS web site, 
these abundance changes should have very small effects on the 
EOS. 



http ://w w w-phy s . llnl . gov/Research/OPAL 
Adopting the usual scale: logn(Li)=12+log(A'Li/A'H)- 
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To determine the extension of the convective envelope the 
classical Schwarzschild criterion has been adopted and the ef- 
ficiency of the energy transport has been modeled with mixing 
length formalism. We did not include rotation. 

Element diffusion is included (Ciacio et al. 1997 with dif- 
fusion coefficients from Thoul, Bahcall & Loeb 1199*31 . while 
radiative acceleration has been neglected since we analysed 
only the PMS phase of low mass stars (M <1.2 M G ). In fact, as 
discussed by Turcotte, Richer & Michaud ( 1998a I and Turcotte 
et al. (1998b), for stars with solar composition and masses 
around 1 M the effect of the radiative levitation becomes com- 
parable to that of the gravitational settling only at the end of the 
MS phase. 

In the following calculations, due to the uncertainties al- 
ready discussed, it is not possible to perform precise evalua- 
tions of the absolute Li surface values; we will instead perform 
a "differential" analysis with respect to a "standard" model, 
investigating the effects of a mixture change. Our "standard" 
model adopts [Fe/H]=0 with the As04 mixture (corresponding 
to a metallicity of Z ~0.013), 7=0.27 and the mixing length 
parameter (a) equal to 1.9. The effects of a variation of [Fe/H], 
Y and a within a plausible range of uncertainty will also be 
discussed. 



3. Li abundance for models with different mixtures 

3.1. Reference models 

We investigated three models of different masses (0.8, 1.0 and 
1 .2 M ), representing the typical mass range covered by most 
observational data for Li in open clusters; higher masses have 
very shallow convective envelopes and thus the effect we are 
studying is too small to be accurately investigated. 

The various models discussed in the present paper are sum- 
marized in Tabled Cols. 1 and 2 report the sequence number 
of the model and the label used in the figures, respectively; in 
Col. 3 the adopted solar mixture is given, while Col. 4 lists the 
variations of elements with respect to the solar mixture - if any 
- or the variations in the surface opacity or global metallicity. 

The upper panel of Fig.^shows the evolutionary tracks for 
the three selected masses from the PMS to the bottom of the 
red giant branch (RGB) with our reference choice of metallic- 
ity and heavy element mixture (model no. 1, As04, in Table 1). 
Note that for completeness we describe here the evolution up to 
the RGB; however, the following discussions will focus only on 
the evolution up to the ZAMS phase where the largest amount 
of Li depletion occurs. As well known, the treatment of the 
early PMS phases is very problematic and the definition of the 
birth line requires the adoption of hydrodynamical codes (see 
e.g. Stahler & Palla 2004 for an extensive analysis). However, 
this is not a crucial point for our purposes since the stellar char- 
acteristic after the birth line are largely insensitive to the pre- 
vious evolutionary history. Moreover, for the selected masses 
the equilibrium is reached very early (ages < 1 Myr), when 
they have not yet started surface Li depletion: indeed, at the 
beginning of the PMS phase, the zone where the Li burning 
temperature is exceeded represents a small fraction of the total 
external convective region; the dilution is thus very large. For 



these reasons we will show our calculations starting from an 
age of 1 Myr. 

The middle panel of Fig.^shows the behavior of Li abun- 
dance {Nu/Nh) as a function of the age for the same models, 
while in the lower panel the temperature at the bottom of the 
convective zone {T cz ) is reported. As well known, the smaller 
the mass the deeper the convective zone and thus the higher the 
temperature at its lower edge. During the PMS the extension of 
the convective envelope and T cz decrease until, in the cases of 
1 and 1.2 M , it becomes lower than ~ 2.5 x 10 6 K; as a con- 
sequence Li destruction stops and the surface abundance de- 
creases only due to microscopic diffusion. On the contrary, for 
the model with 0.8 M Q Li burning is efficient also during the 
MS phase even if the largest amount of Li is depleted during 
the PMS phase (note that MS Li depletion is not evident from 
the figure due to the large scale adopted: for instance, from the 
ZAMS to the turn off, Nu/Nu decreases from ~ 6 x 10~ 12 to 
~ 1.5 X 10~ 13 , i.e. by a factor -40). 

After central hydrogen exhaustion, the H burning continues 
in a shell surrounding the He core and the stars move to lower 
effective temperatures. During this phase the convective enve- 
lope extends into regions in which Li has accumulated due to 
the effect of diffusion, and the external abundance temporarily 
increases; then, it decreases below the level of detection when 
the external convective region reaches the layers in which Li 
has been completely burned (dilution effect). 

The details of Li depletion depend on the adopted physical 
inputs, in particular the radiative opacity and chemical compo- 
sition for a given mixing scheme. Thus, the situation shown in 
Fig- H has to be considered as a description of the dependence 
of Li depletion on the mass, while the detailed quantitative val- 
ues hold only for the labeled parameters. 

3.2. Global uncertainties 

Empirical data and theoretical predictions are usually com- 
pared in the [logn(Li), r e ff] plane; thus it is important to eval- 
uate the ability of canonical models to reproduce these quanti- 
ties. The present empirical uncertainties on the metallicity, He 
abundance and efficiency of the external convection (usually 
parameterized by the mixing length value a) play an important 
role in the comparison between theory and observations. 

At the typical effective temperatures of the stars under in- 
vestigation the surface He abundance cannot be measured spec- 
troscopically. As a consequence one has to rely on indirect 
methods, such as assuming a relation between metal and He 
enrichment of the Galactic medium which is not well defined 



(see, e.g., Pagel & Portinari 1998 , Castellani, DeglTnnocenti & 



Marconi 1999 , Luridiana et al. 2003 , Olive & Skillman 20041. 
Under this assumption, AF=+0.02 is generally assumed as an 
estimate of the original He abundance uncertainty for these 
stars. 

Regarding the Fe content, disk stars cover a wide interval of 
[Fe/H] and typical uncertainties depend on several factors such 
as on whether the analysis is spectroscopic or photometric, the 
quality of the data, the method of analysis, statistics etc. In the 
case of young nearby stars (see e.g. Chen et al. 2000; Bensby 
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Table 1. Summary of the models computed in this work. In all cases the total metallicity is fixed to Z = 0.013, with the exception 
of model no. 16; the PMS evolution has been followed for the three masses 1.2, 1.0 and 0.8 M . 



No. 



Label 



Mixture 



Element, opacity, or metallicity variations 



1 


As04 (standard) 




Asplund et al. (2004) 


2 


GN93 




Grevesse& Noels (1993) 


3 


GS98 




Grevesse&Sauval(1998) 


4 


As05 




Asplund, Grevesse & Sauval 2005 


5 


As04 (atmospheric opacity deer. 


by 5 %) 


Asplund et al. (2004) 


6 


As04 (atmospheric opacity deer. 


by 10 %) 


Asplund et al. (2004) 


7 


OnlyC 




- 


8 


OnlyO 




- 


9 


Only Ne 




- 


10 


C(-20%) 0(+) 




Asplund et al. (2004) 


11 


O(-20%) C(+) 




Asplund et al. (2004) 


12 


Ne(+30%) O(-) 




Asplund et al. (2004) 


13 


Ne(-30%) 0(+) 




Asplund et al. (2004) 


14 


Heavy elements (+20%; 


i 


Asplund et al. (2004) 


15 


CNO (+20%) 




Asplund et al. (2004) 


16 


GN93Z = 0.018 




Grevesse & Noels (1993) 



Standard solar mixture 

Standard solar mixture 

Standard solar mixture 

Standard solar mixture 

Atmospheric opacity decreased by 5 % 

Atm. opac. decreased by 10 % 

Metallicity composed by 100% C 

100% O 

100% Ne 

C mass fraction decreased by 20% (compensated by O enhancement) 

O mass fraction decreased by 20% (compensated by C enhancement) 

Ne mass fraction enhanced by 30% (compensated by O decrease) 

Ne mass fraction decreased by 30% (compensated by O enhancement) 

Elements heavier than O increased by 20% 

CNO increased by 20% 

Z - 0.018, standard solar mixture 



et al. l2003t and open clusters with nearly solar metallicity (e.g. 
Sestito et al. l2003l Carretta et al. l2004l Randich et al. 2006 1 the 
quoted errors range from A[Fe/H]~+0.01 up to ~+0. 06-^0. 07, 
when the analysis is based on a large sample of stars observed 
with high resolution spectrographs. Errors can be as high as 
A[Fe/H]~+0.1 (or more) if results from low resolution spec- 
troscopy or from the photometry are considered (see e.g. Friel 
et al. 120021 Nordstrom et al. 2004). We adopt here a value of 
A[Fe/H]=+0.1 as a conservative estimate of observational un- 
certainties on the metallicity of solar- type stars. 

Finally, the a parameter adopted in the calculations, de- 
pending on the selected set of atmospheric models, could be 
calibrated by requiring the reproduction of the color of the MS 
stars of the observed open cluster. However a is a useful param- 
eter to take into account our lack of precise knowledge about 
the convection mechanisms, and there is no reason for which 
it must assume the same value for stars of different mass and 
chemical composition, or even for different evolutionary phases 
of the same star. When a direct calibration of the a value for 
PMS stars through the reproduction of their stellar colors is not 
possible, one must take into account an uncertainty on this pa- 
rameter; we allow a variation of the mixing length parameter a 
by +0.3. 

We computed six models for each of the three selected 
masses by fixing, in turn, one of the three parameters at its 
boundary values (lower and upper) and keeping the others at 
their central values (Z=0.013, 7=0.27, a— 1.9). In this way, we 
obtained relative variations of Li abundance and effective tem- 
perature with respect to our "standard" model, representing an 
evaluation of the uncertainties on Li and T e g due to plausible 
empirical errors on Z, Y and a. Then, the relative variations of 
Li and T e g (on the increase and on the decrease) due to each of 
the quoted parameters are added quadratically to give a quanti- 
tative estimate of the total (upper and lower) uncertainties. 

In Fig.|5]the results for our "standard" models for the three 
different masses at 30 Myr (when for all the stars the main Li 
depletion has ended) are shown in the [logn(Li), T e ff] plane 
with the error bars due to the ranges of [Fe/H], Y and a, calcu- 
lated as described above. The error bars represent the range of 
predicted theoretical values due to the uncertainties discussed 
above, while intrinsic errors of the model (due to physical in- 



puts) are not considered; therefore the error bars are an indica- 
tive estimate of the uncertainties in the initial physical parame- 
ters. Note that among the three quoted error sources, the largest 
contribution to the total error bar is due to the uncertainty on 
Z; in particular, the effect of the metallicity variation is very 
strong for the lowest mass, for which log «(Li) appears to be 
unpredictable. 

3.3. Heavy element mixture variations and opacity 
effect 

The next step is aimed to evaluate the effect of a variation of 
the adopted heavy element mixture with respect to the sources 
of uncertainty described above. Thus we calculated models for 
the same (Z, Y, a) and the Grevesse & Noels (1993), Grevesse 
& Sauval (1998) and Asplund et al. (2005) mixtures (models 
no. 2, 3, 4 in Table 1; GN93, GS98, As05). Since the temper- 
atures of main interest for Li depletion are those at the bottom 
of the external convective region during PMS Li burning, that 
is T cz ~ 3 -f 4 x 10 6 K (see Fig.0, we calculated high tempera- 
ture opacities (T> 12000 K) for the chosen mixtures using the 
procedures available on-line from the OPAL group. 

The results of our calculations are reported in Fig. 13 The 
variations in Li abundance and 7^ due to the mixture changes 
are within the (conservative) previously quoted uncertainties 
but they are not negligible. It can also be noticed that the 
smaller the mass, the more sensitive is the dependence of the 
Li abundance on the adopted physical parameters. 

For completeness Fig. [3] shows the surface Li abundance 
behavior during the PMS phase for the three selected masses 
with the four different mixtures of Fig.|2] The models that un- 
dergo the largest amount of Li depletion have higher opacity in 
the region at the bottom of the convective zone, and thus deeper 
and hotter bottom convective envelopes (see Sect.^Ji. 

However, the various quoted mixtures differ in the abun- 
dance of many elements. In the following (Sect. |4}, we will 
thus discuss simpler situations in which only one or a few crit- 
ical element abundances are changed. This provides a better 
understanding of the dependence of the opacity at the bottom 
of the convective envelope - and thus of Li depletion - on the 
element abundance variations. 
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Low temperature (i.e. atmospheric) opacities calculated for 
the As04 and As05 mixtures are not available; however, since 
the temperature of the bottom of the convective zone is rela- 
tively high, the effect of a change in the surface opacity due to 
the mixture update on Li depletion should not be too relevant. 
We expect that changes in atmospheric opacities mostly affect 
the model with the smallest mass (0.8 M ), which has the low- 
est effective temperature and the deepest convective envelope. 
From the calculations by the OPAL group we found that, for a 
0.8 M star and at a temperature of about 12000 K (the lowest 
T for which the OPAL tables are computed), the change in the 
internal opacity values when going from the GN93 to the As04 
mixture (at fixed Z) is less than 5%; we took this value as rep- 
resentative also of a possible variation in the lower temperature 
opacities. In order to analyse the effect of such a variation, we 
evolved a 0.8 M (Z=0.013, 7=0.27, a=1.9) with the usual at- 
mospheric opacities (Alexander & Ferguson 1994) artificially 
decreased by 5 and 10% (models no. 5 and 6 in Table 1). For 
the interior opacities we adopted the As04 mixture. The results 
for Li depletion are shown in Fig.0]compared with the "stan- 
dard" 0.8 M Q . The effect is rather small; the decreased atmo- 
spheric opacities result in models with hotter external regions, 
thus, even if the depth of the convective envelope is nearly 
unchanged, their bottom temperature is slightly higher and a 
slightly larger amount of Li depletion occurs. Nevertheless, the 
effect is negligible with respect to that due to a mixture varia- 
tion for the interior: the ZAMS Li abundances (age ~ 30 Myr) 
decrease by factors ~1.5 and 3 (with respect to the As04 stan- 
dard model) when surface opacities are decreased by 5 and 10 
%, respectively; on the other hand, going from the standard 
As04 mixture to the GN93 one, the ZAMS Li content increases 
by a factor ~18. 

4. Interpretation of the Li abundance behavior 

In order to understand the effect of the different mixtures on 
the extension of the convective envelope, we analyse the con- 
tribution of the different elements to the radiative opacity at the 
temperatures of interest, i.e. those at the base of the convective 
envelope during the PMS. Therefore, we computed opacity ta- 
bles for pure elements by using the data and routines available 
on-line form the LAOL group 3 (we did not use OPAL opacities 
in this case because they are not available for a single element). 

Figure [5] shows the pure element opacity for the most im- 
portant elements: H, He, C, O, Ne, Fe, Si. The opacities are 
calculated for T = 3.5 ■ 10 6 K and T = 4.6 • 10 6 K and for 
densities ranging from to 10 g cm -3 . Note that at the selected 
temperatures the most opaque element is Fe, followed by Ne, 
and O is more opaque than C; the opacity of pure He and H is 
negligible with respect to those of metals. Obviously in a real- 
istic case these single element opacities must be weighted with 
the abundance of each element in a mixture in which H and He 
are present too. 

The next step to understand the results of Sect.|3]is mus t0 
compute PMS models in which, for the standard metal and He 
abundances, the contribution to the metallicity is made by only 



http://www.t4.lanl.gov/cgi-bin/opacity/tops.pl 



one element. Therefore, we evolved the models for Z=0.013 
and T=0.27 by assuming that the contribution to the metallicity 
is made by only C, O or Ne; we also took into account the varia- 
tion of the composition in calculating the nuclear reaction rates, 
while for the other thermodynamical quantities we adopted the 
OPAL EOS_2001 tables calculated for the Grevesse (1991) so- 
lar mixture. Our assumptions, as a first approximation, affect 
the PMS evolution only through the variation of the stellar 
opacity; we checked that the temperature profile in the con- 
vective envelope does not change in a significant way for the 
three models. 

The opacity tables were computed by adopting the proce- 
dures available on-line from the OPAL group. The results of 
these numerical experiments are shown in Fig.[6](models no. 7, 
8, 9 of Table 1); the behavior of Li abundance is easily ex- 
plained by the opacity values seen in Fig. [5] the element with 
the highest opacity at temperatures and densities typical of the 
base of the convective envelope leads to the highest Li deple- 
tion for the models of Fig. [6] In addition, the results for our 
"standard" models are not so different from those obtained for 
a metallicity contributed by only O: this is easily understand- 
able, since O is the most abundant element (after H and He) 
and its opacity is greater than that of C (which is the next most 
abundant element). The behavior is similar for the three masses 
although, as already noted in Fig. [3] the effect is smaller for 
larger stellar masses. 

A more realistic situation is represented by models with our 
"standard" (As04) mixture in which only few elements are var- 
ied each time by ~ 20 -f 30%, keeping the total metallicity con- 
stant. This could reflect not only an update of the solar mixture, 
but also a possible scatter in the composition of an open cluster: 
indeed, although stars in the same cluster are usually assumed 
to have the same chemical composition, an intrinsic dispersion 
in [Fe/H] (up to ~ ±0.15 dex) among cluster members cannot 
be excluded (e.g. Nissen 1988; Schuleret al. 2006). Analogous 
variations may exist for other elements. In Fig.0we show only 
the results for 1.0 M Q , since the overall behavior was the same 
as in the previous calculations. The upper panel of Fig.0shows 
the effects on Li depletion when C and O are changed maintain- 
ing the sum unaltered in mass of the two elements (Xc+Xo). In 
the first case C is decreased by 20% in mass while in the sec- 
ond case O is decreased by the same amount (models no. 10 
and 1 1 in Table 1). An enhancement in O abundance increases 
the Li depletion, due to the higher opacity at the bottom of the 
convective zone and thus to the deeper extension of the con- 
vective envelope; the opposite happens when O is decreased. 
However, the effect is not symmetric: this is due to the fact 
that we keep the sum of the abundance in mass of the two el- 
ements - and thus the total metallicity - constant, and the two 
elements have different atomic weights and "standard" abun- 
dances; therefore, C and O are varied by different total amounts 
in number. These variations have different effects on the global 
mixture and on the opacity, and thus on Li depletion. The mid- 
dle panel of Fig.^shows the results of a decrease/increase by 
30% in mass of the Ne mass fraction (with a corresponding O 
variation to preserve Xo+X^ e ; models no. 12 and 13in Table 1). 
Due to the higher Ne opacity, an increase/decrease of Ne abun- 
dance leads to a decrease/increase of the Li surface abundance 



P. Sestito et al.: Heavy element opacity and Li depletion 



but, in this case, the effect is obviously symmetric. The last 
numerical experiment (bottom panel of Fig. [7) evaluates the 
effects of increasing/decreasing the relative abundance of CNO 
elements and elements heavier than O (models no. 14 and 15 
of Table 1): an increase in CNO with respect to heavier ele- 
ments produces a lower depletion, and vice versa. The results 
of all the previous numerical experiments can be understood in 
terms of the radiative opacity of the heavy elements. In partic- 
ular, from the point of view of Li depletion of low mass stars, 
the maximum depth reached by the convective envelope is very 
important, which in turn is sensitive to the opacity in the tem- 
perature range 3-r4- 10 6 K (see Fig.Q. In this interval the most 
opaque elements, among those with a significant abundance, 
are the elements heavier than O. If the abundance of these ele- 
ments is higher, then the mixture is more opaque and thus the 
convective envelope is deeper. At fixed metallicity (Z=0.013), 
the total abundance of the elements heavier than O increases 
progressively in the following mixtures: GN93, GS98, As05 
and As04, thus explaining the behavior shown in Fig. [3] 

5. Models with fixed [Fe/H] 

In the previous sections we analysed the effects of a mixture 
variation at fixed global metallicity: however, the variation in 
the solar composition also affects the conversion of the [Fe/H] 
value of the observed stars (derived through a spectroscopic 
analysis) to the total metallicity Z. This point is very subtle 
because to compute stellar models one needs both the Fe abun- 
dance and the total metallicity Z; the two quantities are related 
through the distribution of elemental abundances. In general, 
unless individual abundances are derived, a solar mixture is 
adopted. As a consequence, a revision of the photospheric so- 
lar abundance directly leads to a variation of the inferred to- 
tal metallicity from the observed [Fe/H]. For solar-like stars, 
[Fe/H] can be estimated by means of differential analysis, that 
is by comparison of the Fe lines of the star with those of the 
Sun. In such a case, the numerical value of [Fe/H] is unaffected 
by a change of the solar mixture and the updated global metal- 
licity of the star can be easily derived by adopting the new value 
of (Z/X) . That is: 

log(Z/X). = [Fe/H] + log(Z/X) . 

In order to give a quantitative estimate of the change of Z 
(at fixed [Fe/H]) when the mixture and Z/X Q are altered, for the 
GN93 composition [Fe/H]=0 corresponds to Z~0.018, while 
for the As04 mixture, as already discussed, the correspond- 
ing metallicity is Z~0.013 (standard value). Figure [8] shows 
the PMS Li depletion for a 1 M model (F=0.27, a=1.9) 
for the GN93 composition and Z=0.018 (model no. 16 in 
Table 1) compared to our "standard" model (As04 mixture and 
Z=0.013); for reference, we also plot the model with GN93 
mixture and Z=0.013, already reported in Fig. [3] For the model 
with Z=0.018 the effect of an increase in the opacity - due 
to the increase of Z - overcomes the influence of the mixture 
change. 

We show in Fig. [5] a comparison between our "standard" 
model (As04 mixture), some of the models presented in Fig. [7] 



(all with Z = 0.013), the GN93 model with Z = 0.018, and the 
observational data for three young open clusters (ages ~30-^50 
Myr), in which solar-type stars are on the ZAMS. All the clus- 
ters plotted have [Fe/H]~0: IC 2602 (Randich et al. I207HI . 
IC 23 91 (Stauffer et al. [T989l, a nd IC 4665 (Martin & Monies 
119971 see also Shen et al. 2005 for the estimate of the metallic- 
ity). The Li data for the clusters have been analysed by Sestito 



& Randich (2005), starting from the equivalent widths pub- 
lished in the literature and using the same method of analy- 
sis. This plot compares the variations of our theoretical results 
- due to changes in the composition - and Li abundances ob- 
served in stars, and is not an attempt to fit the empirical data. 
For the reasons discussed in Sect. 1, even an apparent agree- 
ment between the theoretical predictions and the data of a clus- 
ter with a given age could be fictitious; to claim a real agree- 
ment between theory and observations, i.e. the goodness of the 
physical inputs and the macroscopic mechanisms adopted in 
the calculations, the models should be simultaneously com- 
pared to data from clusters with different chemical composi- 
tions and ages. The plot clearly shows how the match between 
theory and observations is strongly dependent on the adopted 
metallicity mixture and how the situation gets worse towards 
lower masses. 

6. Summary 

We presented an investigation of the effects of the chemical 
abundances on stellar opacity and 7 Li depletion. Within the 
"standard" scenario for stellar evolution - in which the po- 
sition of the convective boundaries are determined following 
the Schwarzschild criterion - the theoretical predictions for the 
temporal behavior of surface Li depletion are affected by vari- 
ous uncertainties; among these a major role is played by radia- 
tive opacity and therefore by the metallicity and the distribution 
of the single elements in the global mixture. 

Focusing on pre-main sequence evolution, where the 
largest amount of Li burning is predicted by standard theories, 
we computed stellar models for three selected masses (0.8, 1.0 
and 1.2 Mq, with Z = 0.013, T=0.27, ar=1.9) by varying the 
global mixture (Grevesse & Noels 1993; Grevesse & Sauval 
1998; Asplund et al. 2004; Asplund, Grevesse & Sauval 2005), 
showing that the effects of the chemical abundances must be 
taken into account. We have quantitatively computed the con- 
tribution of elements heavier than O for Pop. i models. The 
comparison between models with different mixtures revealed 
that, at fixed Z, a larger amount of Li is depleted in the pres- 
ence of a larger fraction of elements heavier than O. Therefore 
these species - although less abundant than CNO - play a fun- 
damental role in stellar opacity and PMS Li evolution. 

These results have been interpreted by analysing the contri- 
bution of the single elements to the opacity at the temperatures 
and densities typical of the base of the convective envelope and 
the effects on Li depletion of single element variations in the 
mixture. 

In addition, we computed models with the GN93 solar com- 
position in which the global metallicity Z is varied, show- 
ing that the effect of increasing Z overcomes the influence of 
the mixture change. We conclude that changes in the chemi- 
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cal composition do affect Li depletion. Comparisons between 
model predictions and observations must take into account the 
element distribution in the stars. Moreover, the abundance de- 
terminations for open clusters should rely on the 3-D atmo- 
spheric models which are more consistent than the usually 
adopted 1-D models. None of the adopted solar mixtures and 
variations in the chemical composition appear to explain the 
very small amount of PMS Li depletion observed for solar-type 
stars. 

Acknowledgements. We are extremely grateful to V. Castellani and S. 
Shore for a careful reading of the manuscript. We warmly thank the 
referee (A. Weiss) for the useful comments. Financial support for this 
work was provided by the Ministero dell'Istruzione, dell'Universita 
e della Ricerca (MIUR) under the scientific project "Continuity and 
discontinuity in the Galaxy formation" (PI.: R. Gratton). 



References 

Alexander, D.R., & Ferguson, J.W. 1994, ApJ, 437, 879 

Asplund, M., Grevesse, N, Sauval, A.J., Allende Prieto, C, & 

Kiselman, D. 2004, A&A, 417, 751 
Asplund, M., Grevesse, N., & Sauval, A.J. 2005, in 

"Cosmic Abundances as Records of Stellar Evolution and 

Nucleosynthesis", eds. FN. Bash, & T.J. Barnes, ASP Conf. 

Series, 336, 25 
Bahcall, J.N., & Serenelli, A.M. 2005, ApJ, 626, 530 
Bensby, T., Feltzing, S., & Lundstrom, I. 2003, A&A, 410, 527 
Cariulo, P., DeglTnnocenti, S., & Castellani, V. 2004, A&A, 424, 927 
Carretta, E., Bragaglia, A., Gratton, R., & Tosi, M. 2004, A&A, 422, 

951 
Castellani, V., DeglTnnocenti, S., & Marconi, M. 1999, MNRAS, 303, 

265 
Chen, Y.Q., Nissen, P.E., Zhao, G., Zhang, H.W., & Benoni, T 2000, 

A&AS, 141,491 
Chieffi, A., & Straniero, O. 1989, ApJS, 71, 47 
Ciacio, F, DeglTnnocenti, S., & Ricci, B. 1997, A&AS, 123, 449 
D' Antona, F, & Mazzitelli, I. 1997, MmSAIt, 68, 807 
Friel, E., Janes, K.A., Tavarez, M., et al. 2002, AJ, 124, 2693 
Grevesse, N. 1991, A&A 242, 488 
Grevesse, N., & Anders, E. 1991, in "Solar Interior and Atmosphere", 

Tucson, AZ, University of Arizona Press, p. 1227 
Grevesse, N., & Noels, A. 1993 in "Origin and Evolution of the 

Elements", eds. N. Prantzos, E. Vangioni-Flam, & M. Casse, 

Cambrdige University Press, Cambridge, p. 15 
Grevesse, N, & Sauval, A.J. 1998, SSRv, 85, 161 
Iglesias, C, & Rogers, F.J. 1996, ApJ, 464, 943 
Jeffries, R.D. 2000 in "Stellar Clusters and Associations: Convection, 

Rotation, and Dynamos", eds. R. Pallavicini, G. Micela, & S. 

Sciortino, ASP Conf. Ser. 198, p. 245 
Jeffries, R.D. 2006, in "Chemical Abundances and Mixing in Stars in 

the Milky Way and its Satellites", eds. L. Pasquini, & S. Randich, 

ESO Astrophysic Symposia, in press 
Luridiana, V., Peimbert, A., Peimbert, M., & Cervino, M. 2003, ApJ, 

592, 846 
Martin, E., & Montes, D. 1997, A&A, 318, 805 
Nissen, P. 1988, A&A, 199, 146 

Nordstrom, B., Mayor, M., Andersen, J., et al. 2004, A&A, 418, 989 
Olive, K.A., & Skillman, E.D. 2004, ApJ, 617, 29 
Pagel, B.E.J., & Portinari, L. 1998, MNRAS, 298, 747 



Palla, F 2002 in "Star Formation and the Physics of Young Stars", 

eds. Bouvier, & J. -P. Zahn, EDP Sciences, Aussois (France) EAS 

Publications Series, Vol. 3 p. Ill 
Piau, L. 2005, submitted to ApJ, astro-ph/05 1 1402 
Piau, L., & Turck-Chieze, S. 2002, ApJ, 566, 419 
Piau, L., Randich, S., & Palla, F 2003, A&A, 408, 1037 
Pinsonneault, M. 1997, ARAA, 35, 557 
Potekhin, A. Y. 1999, A&A, 351, 787 
Randich, S. 2006, in "Chemical Abundances and Mixing in Stars in 

the Milky Way and its Satellites", eds. L. Pasquini, & S. Randich, 

ESO Astrophysic Symposia, in press 
Randich, S., Pallavicini, R., Meola, G., Stauffer, J.R., & Balachandran, 

S.C. 2001, A&A, 372, 862 
Randich, S., Sestito, P., Primas, F, Pasquini, L., & Pallavicini, R. 

2006, in press, astro-ph/0601239 
Rogers, F.J., 2001, Contributions to plasma physics, V41(N2-3), 179 
Schuler, S.C, Hatzes, A..P, King, J.R., Kuerster, M., & The, L.-S. 

2006, AJ, 131, 1057 
Sestito, P., & Randich, S. 2005, A&A, 442, 615 
Sestito, P., Randich, S., Mermilliod, J.-C, & Pallavicini, R. 2003, 

A&A, 407, 289 
Shen Z.-X., Jones, B., Lin, D.N.C., Liu, X.-W., & Li, S.-L. 2005, ApJ, 

635, 608 
Siess, L., Dufour, M., & Forestini, M. 2000, A&A, 358, 593 
Stahler, S.W., & Palla, F 2004, "The Formation of Stars", New York, 

Wiley 
Stauffer, J., Hartmann, L.W., Jones, B.F, & McNamara, B.R. 1989, 

ApJ, 342, 285 
Swenson, F.J, Faulkner, J, Iglesias, C.A., Rogers, F.J., & Alexander, 

D.R. 1994a, ApJ, 422, L79 
Swenson, F.J, Faulkner, J, Rogers, F.J., & Iglesias, CA. 1994b, ApJ, 

425, 286 
Thoul, A.A., Bahcall, J.N., & Loeb, A. 1994, ApJ, 421, 828 
Turcotte, S., & Christensen-Dalsgaard, J. 1998, Space Sci. Rev., 85, 

133 
Turcotte, S., Richer, J., & Michaud, G. 1998a, ApJ, 504, 559 
Turcotte, S., Richer, J., Michaud, G., Iglesias, C, & Rogers, F 1998b, 

ApJ, 504, 539 



Z=0.013Y=0.27 




P. Sestito et al.: Heavy element opacity and Li depletion 

4 



3 - 



Asplund et al. 2004 (As04) mixture 



3.75 3.7 

logfUK]) 




Z=0.01 3 Y=0.27 As04 mixture 

1.2M, 

TO 



0.8M,„ 




log (age [yr]) 





1.0M„ 



c 2 



1 - 



• As04 
■ As05 
^GS98 
*GN93 



0.8M = , 



6800 6400 6000 5600 5200 4800 4400 4000 
T et ,[K] 

Fig. 2. Theoretical Li abundance (expressed as logn(Li) = 
12 + log(Nu/Nn)) as a function of the effective temperature 
for 0.8, 1.0 and 1.2 M at the age of -30 Myr. Results for the 
four labeled different mixtures are shown (Grevesse & Noels 
1993 - GN93; Grevesse & Sauval 1998 - GS98; Asplund et 
al. 2004 - As04, our "standard" mixture; Asplund et al. 2005 
- As05). For As04 we also show the error bars due to the esti- 
mated uncertainties on [Fe/H], Y and a (see text). 



log (age [yr]) 

Fig. 1. Upper panel: the evolution in the HR diagram of our 
three selected "standard" models (0.8, 1.0 and 1.2 M Q , with 
Z=0.013, 7=0.27 and or=1.9; Asplund et al.2004 - As04 - 
mixture) from the PMS to the RGB phase. Middle panel: sur- 
face Li abundance as a function of time (log age) for the three 
models. The abundance is expressed as the ratio between Li and 
H numerical abundances (Nu/Nu)', the starting of the MS phase 
and the turn off (TO) positions are indicated. Lower panel: tem- 
perature at the bottom of the convective envelope (T cz ) vs. log 
age (when the star is fully convective the bottom of the convec- 
tive core coincides with the center of the star); the indicative 
temperature of ignition for the 7 Li+p reaction is ~ 2.5 x 10 6 K. 
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Fig. 3. Surface Li abundance behavior during the PMS phase 
for the different labeled mixtures. The three selected masses 
are reported: 0.8 M (upper panel), 1.0 M Q (middle panel), 1.2 
M Q (lower panel), for Z=0.013, 7=0.27, a=1.9. Abundances 
and ages are expressed as in Fig.^ note the different scales for 
Li abundances adopted in the three panels. 
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Fig. 4. Surface Li abundance vs. age for 0.8 M (Z=0.013, 
7=0.27, a=1.9). The standard As04 model is compared to 
models in which the atmospheric opacity has been decreased 
by 5% (dotted line) and by 10% (long dashed line). We report 
for reference also the model with the mixture by GN93. 




Fig. 5. The Rosseland opacity for T = 3.5 ■ 10 6 K 
(left panel) and T - 4.6 • 10 6 K (right panel) as a 
function of the density for several pure elements, as la- 
beled. The opacity values are obtained from the URL 
http://www.t4.lanl.gov/cgi-bin/opacity/tops.pl 
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Fig. 7. Li abundance vs. time for PMS models of 1 M Q . We 
Fig. 6. Surface Li abundance as a function of time for PMS adopted , as usual, Z=0.013, Y=0.27, and a=l.9 and we var- 
models of 0.8 (upper panel), 1 .0 (middle panel) and 1 .2 (lower ied the mass fraction of some elements starting from me As04 
panel) M , For each mass Li abundance is calculated for mixture (see text) 
Z=0.013, Y-0.21, a-\.9 for a metallicity due to only one ele- 
ment: C, O, and Ne. The results are compared with the "stan- 
dard" model, where the As04 mixture is adopted. Note the dif- 
ferent scales adopted in the three panels. 
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Fig. 8. Surface Li abundance vs. age for our standard model 
(As04 mixture with Z = 0.013) compared to a model with Z = 
0.013 and GN93 mixture (the same of Figs. 00 and© and 
to a model with the GN93 solar composition and metallicity 
Z = 0.018 (see text). 
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Fig. 9. Li abundance vs. T e e: comparison between the obser- 
vational data for young open clusters (ages ~30-h50 Myr) and 
the predictions of models. For the sake of clearness the calcu- 
lations done for the three selected masses for each mixture are 
connected by lines. The mean observational error bars for the 
clusters are also shown. 



